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Using the numerical diagonalization method, we investi- 
gate the one-dimensional d-p model, simulating a Cu-O lin- 
ear chain with strong Coulomb repulsions. Paying attention 
to the effect of the transfer energy t pp between the nearest 
neighbor oxygen-sites, we calculate the critical exponent of 
correlation functions K p based on the Luttinger liquid rela- 
tions and the ground state energy Eo(4>) as a function of an 
external flux <j>. We find that the transfer t pp increases the 
charge susceptibility and the exponent K p in cooperation with 
the repulsion Ud at Cu-site. We also show that anomalous flux 
quantization occurs for K p > 1. The superconducting region 
is presented on a phase diagram of Ud vs. t pp plane. 

PACS numbers: 71.10.Fd, 71.10.Hf, 74.20.Mn 

Low-dimensional strongly correlated electron systems 
attract much interest due to the possible relevance to 
high-T c superconductivity. A considerable number of 
theoretical studies have been performed on the strongly 
correlated fermion models either in one-band or multi- 
band version to understand the mechanism of the high 
temperature superconductivity. However, there is no 
consensus on what the actual mechanism of the super- 
conductivity is, except for the fact that the strong repul- 
sive interaction on Cu-site plays an important role in the 
mechanism. |l],^ To achieve a solid understanding of this 
problem, some rigorous studies, even of the simplest mod- 
els, are highly desirable at this stage. For such purpose, 
numerical diagonalization studies of finite size systems 
may be successful. In particular, one-dimensional 

d-p model, simulating a Cu-0 linear chain with strong 
Coulomb repulsions, provides a good target to investigate 
the mechanism of the high-T c superconductivity. j?] 
The model contains hopping t p d between Cu(d)-site and 
oxygen(p)-site and repulsive interaction at the d- and p- 
site (Ud and U p , respectively). In addition, it can contain 
the nearest-neighbor d-p interaction U p d and/or hopping 
t P p between the nearest-neighbor p-sites. 

In the previous work, the present authors have 

studied the one-dimensional d-p model with large on-site 
Coulomb repulsion Ud at Cu-sites and inter-site repul- 
sion Upd, by using the numerical diagonalization method. 
With the help of the Luttinger liquid relations, the su- 
perconducting(SC) correlation is found to be dominant 
compared with the CDW and SDW correlations in the 
proximity of the phase separation. The SC phase ap- 
pears for 1.0 ^ Upd ~ 1-6 almost independent of the 
filling with 1.3 n < 1.7. Near the half-filling n = 1 
and the full-filling n = 2, however, the SC correlation is 



suppressed. [jl0| Recently, Sudb0 et al. showed that the 
d-p chain with large U p d exhibits flux quantization with 
charge 2e and slow algebraic decay of the singlet SC cor- 
relation function on oxygen sites, jllj Other studies of 
the d-p chain also showed that U p d enhances the charge 
fluctuation as well as the SC fluctuations and when it ex- 
ceeds a certain critical value the system becomes unsta- 
ble towards phase separation. Jl2|,[l3| All of these studies 
claimed that the parameter U p d is central in inducing a 
SC state at T = 0. 

However, if one adds the hopping term t pp to the 
model, the situation will be changed completely. It en- 
hances the charge fluctuation and increases the exponent 
K p as the parameter U p d does. Therefor, the repulsive 
interaction U p d is not always necessary for the SC state. 
In the present paper, to clarify the effect of t pp , we study 
the one-dimensional d-p model with the hopping t pp by 
the numerical diagonalization. We calculate the criti- 
cal exponent of correlation functions K p and the ground 
state energy Eo(4>) as a function of an external flux <f>. 
Using the results for the critical exponent, we show a SC 
region on the Ud — t pp plane. 

We consider the following model Hamiltonian for the 
Cu-0 chain: 

H = t pd (Phdj" + h - c -) + hp (pLpj° + h -c) 

<ij>,<7 <ij>,o 
+ €d H d ]a d 3° +e pJ2 Pi^P™ + Ud J2 n m n dH (1) 

where d^ a and p\ a stand for creation operators of a hole 
with spin a in the Cu(d)-orbital at site j and of a hole 
with spin a in the 0(p)-orbital at site i respectively, and 
ndja = dj a djcr- Here, t p d stands for the transfer energy 
between the nearest neighbor d- and p-sites, which will 
be set to be unity (t p d=V) hereafter in the present study. 
The atomic energy levels of d- and p-orbitals are given 
by td and e p , respectively. The charge-transfer energy A 
is defined as A = e p — e<j- 

To achieve systematic calculation, we use the periodic 
boundary condition for Nh = 4m + 2 and antiperiodic 
boundary condition for Nh — 4m, where Nh is the to- 
tal hole number and m is an integer. This choice of 
the boundary condition removes accidental degeneracy 
so that the ground state might always be a singlet with 
zero momentum. The filling n is defined by n — Nh/N u , 
where N u is the total number of unit cells ( each unit 
cell contains a d- and a p-orbitals) and the Fermi wave 
number kp is given as kp = \n. 



1 



We numerically diagonalize the Hamiltonian up to 12 
sites (6unit cells) using the standard Lanczos algorithm. 
The chemical potential ^{Nh, N u ) is defined by 

(AT AT , E (N h + l,N u )-E (N h -l,N u ) 

where Eo(Nh, N u ) is the total ground state energy of a 
system with N u unit cells and Nh holes. When the charge 
gap vanishes in the thermodynamic limit, the uniform 
charge susceptibility Xc is obtained from 

Xc{N h ,N u ) = — 2 / Nu — — (3) 

fi{N h + 1, N u ) - fj,{N h - 1, N u ) 

In the Luttingcr liquid theory, some relations have 
been established as universal relations in the one- 
dimensional models. p^-[l7|] In the model which is 
isotropic in spin space, the critical exponents of various 
types of correlation functions are determined by a single 
parameter K p , which is the exponent of the power-law 
decay of correlation functions. For the one-dimensional 
d-p model in the weak coupling regime, Matsunami and 
Kimura |^(| calculated the critical exponents of some 
correlation functions by using the renormalization group 
analysis (<?-ology), and shown that the system belongs 
to the Tomonaga-Luttinger (TL) liquid. In the present 
authors' recent work, Q we have shown the numerical re- 
sults of the critical exponent K p and the phase diagram 
on the Ud-Upd plane including the weak coupling region 
and the strong coupling region. In the weak coupling re- 
gion, it agrees with the results of the g-ology. An analysis 
of spin-gap indicates that the strong coupling region also 
belongs to the TL regime and not to the Luther-Emery 
(LE) regime. Here, the LE regime is characterized by 
a gap of the spin excitation spectrum, while in the TL 
regime, the excitation is gapless. 

It is predicted that SC correlation is dominant for 



K p > 1 (the correlation function decays as ■ 



the TL regime and as ~ r K ? in the LE regime), whereas 
Charge Density Wave (CDW) or Spin Density Wave ( 
SDW) correlations is dominant for K p < 1 (correlation 
functions decay as ~ r~( 1+Kp ' in the TL regime and as 
~ r~ Kp in the LE regime). In the TL regime, the singlet 
SC and triplet SC correlation functions have the same 
critical exponent apart from a logarithmic correction. On 
the other hand, in the LE regime, the correlation func- 
tions of SDW and triplet SC decrease exponentially. The 
parameter K p is related to the charge susceptibility y e 
and the charge velocity v c by the relations, [Q,[l5 17J19] 



Kn 



; v cXc 



(4) 



with v r = ^(E 



It: 



i — Eq) , where E\ — Eq is the lowest 
charge excitation energy. We can also determine the K p 
by the Drude weight D 



with D = jj- d , where E(<f) is the total energy of 
the ground state as a function of flux <fi. Jl^] Using these 
two independent equations of the K pi we will check the 
consistency of the Luttinger liquid relations. 

For a non-interacting system (Ud = 0), the Hamilto- 
nian (1) is easily diagonalized and a dispersion relation 
is given as 



E 



= - { €d + e P + 2t pp cos k 



± y / (A + 2t w cosfc) 2 + 16(t pd cos (A;/2)) 2 | (6) 

where k is a wave vector and E + (k)(E~ (k)) represents 
a upper (lower) band energy. Compared with the case 
with t pp = 0, the width of the lower band E~{k) be- 
comes narrower for t pp > 0. Note that it decreases with 
A and becomes perfectly flat at A = 2t pp —t pd /tp P . When 
one decreases A further, the band bends with a peak at 
k = 0. To investigate the band structure of interact- 
ing systems, we consider the chemical potential u, which 
is corresponding to Fermi energy Ep at T = 0. Fig- 
ure ^(a) shows fi as a function of the filling n for sev- 
eral values of Ud with A = 2 and t pp = 0.5. We calcu- 
late the value of /j, by using eq.(2) at the hole densities 
n = 5/4, 7/4, 7/5, 9/5,, 7/6, 9/6 and 11/6. Data points 
with Ud = seem to sit closely along on the dispersion 
curve of E~(kp). It suggests that the finite-size effect is 
sufficiently small. 

When Ud is large, /i is large and weekly depends on n. 
Especially, the slope of /x becomes almost flat at Ud — 6.5. 
In this case, the charge susceptibility Xc almost diverges, 
because of the relation, xj 1 = The result suggests 
that the width of effective hole band becomes very nar- 
row due to a strong correlation effect of large repulsion 
Ud- In order to see the behavior of the charge susceptibil- 
ity more qualitatively, in Fig. |l|(b), we show the inverse 
of the charge susceptibility x as a function of Ud in a 



system with 6 unit cells for A = 2. In the case of t p 



0, 



the [/^-dependence of x _1 is week and x _1 remains as a 
finite value. On the other hand, it is close to zero at 
Ud — 6.75 for tpp = 0.5. It shows that the transfer t pp in- 
creases the charge susceptibility in cooperation with the 
repulsion Ud- If we use the Hartree-Fock(H.F.) approxi- 
mation, the atomic levels of (i-orbitals td are pushed up 
as id — £rf + \Ud < nd >, where < rid > is the mean value 
of hole density at a d-site and determined by solving 
the self-consistent equation of < rid >■ Then, the band 
structure is approximately described as a renormalizcd 
band E{k) whose A is replaced by A = A — \ Ud < nd >■ 
The band width of E~ (k) decreases as A decreases and 
becomes zero at sufficiently large Ud- The decreases of 
the band width leads to the enhancement of Xc as shown 
in Fig. 1(b). 
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FIG. 1. (a) The chemical potential n as a function of the 
filling n for various values of Ud- fl is calculated for N u = 8, 10 
and 12. The solid line represents E~(kp)- The broken lines 
are guide for eye. (b) The inverse charge susceptibility x _1 
as a function of Ud for t vp = 0.0 and 0.5 with the results in 
the H.F. approximation. 

Next, we show the critical exponent of correlation func- 
tions K p as a function of Ud in Fig. 2(a). We calculate 
K p in two ways by using eq.(4) and eq.(5) independently, 
one is obtained by combining the charge velocity v c and 
the charge susceptibility Xc and the other is expressed by 
the Drude weight. The result shows that the values of 
Kp obtained by in the two different ways are well con- 
sistent with each other except for Ud ~ —1- We also 
show K p obtained by the g-ology together in Fig. 2(a). 
In the weak coupling regime, the results of numerical di- 
agonalization are in good agreement with those of the 
g-ology. These results indicate that the Luttinger liquid 
theory can be applied to the one-dimensional d-p model 
successfully and finite-size effect of K p is not so large. 

To judge which of the two estimations of K p is better 
from a point of view of finite-size effect, we watch the 
values of K p at Ud — 0. The one from eq.(4) is equal to 
0.972 and the other from eq.(5) is 1.002. Taking account 
of K p being always equal to unity in noninteracting case, 
we infer that the estimation from eq.(5) is better than 
the other. Furthermore, the comparison of those values 
of Kp with the values obtained in the g-ology in the weak 
coupling regime also indicates that the estimation of K p 
from eq.(5) gives a better result. Therefore, we will use 
the value of K p expressed by eq.(5) in the following. 

Fig. 2 (a) shows that K p decreases as Ud increases when 



Ud ~ 3. However, K p increases with Ud for Ud ~ 3 and 
diverges at Ud ~ 6.8 where the charge susceptibility di- 
verges. The divergence indicates that the width of the 
effective band E~{k) is close to zero. When K p is larger 
than unity, the SC correlation is expected to be most 
dominant compared with CDW and SDW correlations. 
The region where K p > 1 appears at 6.5 ^ Ud ~ 6.8. It 
seems to be fairly small, but surely exists. To confirm the 
ground state to be in the SC state, we show the ground 
state energy E (4>) as a function of an external flux <p 
in Fig. 2(b). The anomalous flux quantization occurs at 
U d = 6.7, where K p is about 2.0. When U d = 6.0, K p 
is less than unity and the anomalous flux quantization 
is not found. These results support that the SC phase 
really appears at K p > 1. 

Generally speaking, the SC phase is induced by an at- 
tractive interaction between quasi-particles. To examine 
the effective attraction, we introduce the repulsive in- 
teractions Udd (U pp ) between the nearest-neighbor d-d 
(p-p) sites and the repulsion U p . If those repulsions act 
between a relevant pair of holes or cancel the effective at- 
traction, the SC state is expected to be suppressed and 
K p is reduced. We find that Udd does not affect K p , but 
Up (> 0.2), as well as U pp , decreases it. It suggests that 
the relevant pair of holes sit not on the d-sites but on the 
p-sites and the strength of the attraction is of order of 
0.2. 
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FIG. 2. (a) K p as a function of Ud- The solid circles 
with the solid line represent the estimation from eq.(5) and 
the open circles with the broken line represent the estimation 
from eq.(4). (b) The energy difference Eo((j>) — -Eo(O) as a 
function of an external flux0. 
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FIG. 3. Contour map for the K p on the Ud-t pp plane at 
A = 2 for the 5-unit 8- hole system (n=8/5). The broken line 
represents t pp obtained in the H.F. approximation. 

Finally, we investigate the critical exponent K p in de- 
tail and obtain the phase diagram in the Ud — t pp pa- 
rameter plane. Figure 3 shows a contour map of K p for 
n = 8/5 system. The region sandwiched between two 
lines with K p = 1 and K p = oo is corresponding to the 
SC region. The region is very small for small Ud, how- 
ever, it increases with Ud- 

In the H.F. approximation, the renormalized band 
E~ (k) becomes flat at a critical value t c vp which is given 
as 

tC P P = \{ A - l U d < > 

+ ^(A-±Ud<n d >r + 8tl d }. 

Here, the approximation breaks down for t pp > tp p . We 
expect that the point of t pp = t^ p corresponds to that 
of K p = oo. In the noninteracting system (Ud = 0), 

tp p = for A = 2. As Ud increases, tp p decreases as 

shown in Fig. 3. It seems that the line of t pp is near that 
of K p = oo. 

In summary, we have numerically diagonalized the one- 
dimensional d-p model with finite sizes. Paying special 
attention to the role of t pp , we have calculated chemical 
potential and critical exponent K p . We have also showed 
the anomalous flux quantization observed in the SC state. 
Using the Luttinger liquid relations, we have found that 
the SC correlation is dominant compared with the CDW 
and SDW correlations in the proximity of the parameter 
region where the renormalized band becomes flat. We 
have confirmed that t pp enhances the charge fluctuation 
and promotes the SC correlation as U p d does. 
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